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A  MODEL  FOR  TARGET  TRACKING  IN  THE 
PRESENCE  OF  GEOGRAPHIC  CONSTRAINTS 


INTRODUCTION 


Incorporating  geographic  constraints  and  environmental  information  into  the  target  tracking 
process  is  a  nontrivial  problem.  To  obtain  accurate  target  estimates,  accurate  signal  modeling 
and  pertinent  state  constraints  must  both  be  present  in  a  fully  integrated  algorithm 
Preconditioning  data  or  applying  state  constraints  after  estimation  does  not  yield  mathematically 
optimal  state  estimates  for  the  typical  problems  encountered. 

Operations  in  littoral  waters  demand  the  ability  to  accurately  track  contacts  in  the  presence 
of  landmass  constraints.  Minefields,  under  appropriate  conditions,  can  also  strongly  bias  the 
probability  of  a  target  presence.  Weather,  sea  state,  and  visibility  can  constrain  or  severely  limit 
target  velocity  and  maneuverability.  Because  the  Navy’s  submarine  focus  has  been  on  open- 
ocean  operations  and  antisubmarine  warfare  (ASW),  these  environment-related  tracking  issues 
have  received  little  attention.  While  some  Navy  tracking  products  such  as  the  Composite  Area 
Analysis  Module  (CAAM),  the  ASW  Tactical  Decision  Aid  (ASWTDA),  and  Visual  Integration 
of  Scene  Information  Over  Networks  (VISION)  allow  for  the  storage,  retrieval,  and  display  of 
environmental  information  (e.g.,  charting,  bathymetry,  sound  velocity  profiling,  shipping  noise 
statistics)  (references  1-3),  current  target  motion  analysis  algorithms  generally  do  not  utilize 
environmental  information.  In  some  instances,  multiple  hypothesis  techniques  are  employed  to 
make  a  propagation  path  decision  when  the  true  path  is  unknown;  yet,  multipath  effects  and  the 
vagaries  of  sound  propagation  in  shallow  water  remain  poorly  integrated  into  the  contact 
management  function. 

As  an  initial  step  toward  fully  incorporating  environmental  information  into  target  motion 
analysis,  this  report  formulates  a  non-homogeneous  process  noise  model  for  application  to  target 
tracking  near  landmasses.  Some  of  the  principal  considerations  in  developing  such  a  model  are 
discussed,  and  the  results  of  scaled  simulations  are  presented  to  illustrate  the  model’s  effect  on 
tracking  algorithm  performance  under  varying  data  conditions. 
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SYSTEM  CONSIDERATIONS 


POINT  ESTIMATE  CONSTRAINTS 

A  fully  coupled,  fully  integrated  target  tracking  algorithm  is  most  desirable,  as  one  is  often 
faced  with  integrating  target  solutions  from  diverse  and  independent  sources.  Even  when 
embedded  in  solutions  presented  by  sensor  data  processing,  kinematic  constraints  induced  by  the 
environment  must  be  considered  when  integrating  solutions  that  may  be  time-late,  aged  by 
prediction,  or  of  uncertain  target  association.  Inclusion  of  environmental  constraint  models 
becomes  even  more  crucial  when  a  predictive  capacity  is  required,  such  as  when  searching  for  a 
lost  contact,  or  when  performing  in-situ  tactical  planning.  Current  tracking  algorithms — if  they 
consider  environmental  constraint  information  at  all — typically  treat  such  information  as  hard 
limits  on  traditional  point  estimates  of  target  state,  with  little  regard  for  the  ramifications  for  the 
state  probability  distribution. 

One  straightforward  method  for  constraining  the  point  estimate  in  an  iterative  algorithm  is 
simply  to  modify  the  magnitude  of  the  state  update  in  each  iteration  to  ensure  constraint 
adherence.  This  approach  is  perhaps  the  least  desirable  means  of  constraint  implementation 
since  the  state  update  prediction  does  not  incorporate  the  constraint,  and  convergence  to  a 
suboptimal  solution  (non-extrema  of  the  constrained  performance  index)  is  possible.  In  addition, 
the  uncertainty  estimate  is  not  affected  by  this  method  of  constraint  implementation;  in  fact, 
when  a  constraint  is  active  (the  estimate  lies  on  the  constraint  boundary),  a  large  portion  of  the 
solution  uncertainty  region  may  violate  the  constraint.  Morevoer,  constraining  point  estimates  in 
this  way  can  lead  to  significant  prediction  errors. 

A  more  rigorous  method  utilizes  penalty  functions  and  Lagrange  multipliers  (reference  4). 
With  this  method  the  constraint  is  incorporated  in  the  state  update  calculation,  and  under 
appropriate  conditions  an  optimal  solution  (an  extrema  of  the  corresponding  cost  function)  is 
ensured.  When  constraint  conditions  are  active,  the  level  curves  of  the  modified  cost  function 
may  provide  qualitative  indications  of  the  solution  uncertainty  region.  However,  since  such 
methods  are  developed  for  functional  minimization/maximization,  they  provide  no  direct  means 
of  quantitative  uncertainty  characterization.  Although  the  cost  function  may  give  an  indication 
of  the  corresponding  uncertainty  region,  there  remains  the  problem  of  characterizing  point 
estimates  near  (but  not  on)  the  constraint  boundary.  Here,  the  inactive  constraints  will  not  be 
incorporated  in  the  cost  function  and  thus  not  properly  accounted  for  in  the  uncertainty  display. 
As  a  result,  parts  of  the  uncertainty  region  may  violate  the  constraint  and,  when  constraints 
become  active,  abrupt  changes  can  occur  in  the  cost  function  and  corresponding  solution 
uncertainty  assessment.  Moreover,  the  additional  processing  required  to  properly  handle  the 
Lagrange  multipliers  at  each  iteration  is  computationally  nontrivial.  Unconstrained  penalty 
function  methods  approximate  the  constrained  functional  minimization/maximization  problem 
by  a  related  unconstrained  problem.  Although  the  discontinuity  at  the  constraint  boundary  is  not 
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a  problem  in  this  case,  the  lack  of  quantitative  uncertainty  measures  and  the  lack  of  appropriate 
target  motion  models  remain  as  issues. 

CONSTRAINED  UNCERTAINTY  REGIONS 

In  the  presence  of  constraints,  full  target  estimate  characterization  requires  the  treatment  of 
non-Gaussian  target  uncertainty  descriptions.  A  constrained  target  motion  model  must  utilize  the 
non-homogeneous  state  transition  probabilities  and  nonstationary  errors  that  characterize  target 
motion  near  constraint  boundaries.  Rigorous  modeling  typically  requires  numerically  intensive 
techniques. 

An  analog  of  the  limited  state  update  is  the  renormalized,  windowed,  probability  density 
function  (PDF)  uncertainty  description.  Here,  the  target  PDF  is  developed  in  the  absence  of  any 
kinematic  constraints  and,  typically,  Gaussian  assumptions/approximations  are  used.  This  PDF 
is  then  multiplied  by  a  window  function  that  is  unity  over  regions  of  navigable  water  and  zero 
where  constraints  are  violated  (i.e.,  over  land).  The  result  must  then  be  renormalized  to  account 
for  the  lost  probability  mass  in  the  inadmissible  region.  Such  a  representation  of  target  state 
probability  will  not  violate  the  constraint  but  may  not  be  representative  of  reality  either, 
especially  when  searching  for  a  lost  target  near  a  peninsula.  The  combination  of  a  point  estimate 
constraint  and  a  windowed  PDF  is  a  crude  but  effective  model  for  the  constrained  state 
estimation  process. 

Because  of  the  nonlinear,  non-Gaussian  nature  of  the  problem,  theoretically  rigorous 
modeling  is  generally  achieved  by  using  numerical  techniques.  Grid-based  methods  maintain 
numerical  values  of  the  target  state  probability  density  in  an  array  representing  a  grid  of  points  in 
the  state  space.  The  equations  for  predicting  a  target  location  or  incorporating  new  measurement 
information  are  not  solved  in  closed  form.  Instead,  the  PDF  is  simply  evaluated  at  each  grid 
vertex  based  on  past  conditional  probability  density  descriptions  evaluated  on  similar  grids.  A 
point  solution  is  extracted  from  the  grid  based  on  the  chosen  metric,  or  the  solution  is 
interpolated  using  the  surrounding  grid  vertices.  Marginalization  is  achieved  through  summation 
over  the  appropriate  indices,  and  uncertainty  regions  are  determined  by  level  curves  (contour 
maps)  of  the  PDF  on  the  grid.  The  difficulties  with  such  methods  are  their  sensitivity  to  grid 
resolution  and  their  computational  burden.  Adaptive  adjustment  of  the  grid  location  and 
resolution  is  helpful  but  carries  a  computational  burden  of  its  own.  If  a  significant  percentage  of 
the  density  lies  between  adjacent  grid  vertices,  or  beyond  grid  boundaries,  information  will  be 
lost,  and  information  lost  due  to  a  poorly  maintained  grid  is  often  difficult  to  retrieve.  An 
adaptive,  grid-based  method  is  employed  in  the  NODESTAR  algorithm  developed  by  Metron 
Inc.  (reference  5). 

Monte  Carlo  methods  avoid  the  need  to  manage  a  grid.  Instead,  the  target  state  PDF  is 
represented  in  histogram  fashion  by  a  large  number  of  points  (“estimates”)  in  the  state  space. 
These  points  may  be  thought  of  as  the  states  of  notional  targets,  and  each  of  these  notional 
targets  evolves  according  to  a  statistical  process  model  that  accounts  for  the  constraint 
(landmass)  presence.  The  probability  distributions  governing  the  choice  of  target  velocity  for 
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propagation  in  time  can  be  adapted  to  include  a  geographic-constraint-avoidance  model.  Such 
methods  are  not  particularly  suited  to  incorporating  new  information  but  are  adept  at  predicting  a 
target  PDF  from  a  given  reference  time  to  some  future  time,  in  light  of  the  given  constraints.  As 
with  grid-based  methods,  there  is  a  significant  computational  burden  associated  with  Monte 
Carlo  techniques  because  of  the  large  number  of  data  points  required  to  generate  a  histogram 
representing  the  target  distribution.  A  Monte  Carlo  tracker  applied  to  the  landmass  avoidance 
problem  is  currently  under  investigation  at  the  Naval  Surface  Warfare  Center  (reference  6). 

Implicit  in  both  the  grid  and  Monte  Carlo  methods  are  statistical  descriptions  of  target 
motion  in  the  proximity  of  the  land  boundary.  Should  the  predicted  position  violate  the  landmass 
constraint  during  the  next  prediction  interval,  remedial  action  is  taken.  In  the  case  of  the  Monte 
Carlo  method  of  reference  6,  notional  targets  that  violate  the  constraint  are  replaced  with  samples 
that  do  not  violate  the  constraint.  The  actual  characterization  of  the  velocity  distribution  from 
which  samples  are  drawn  can  be  tailored  to  the  operating  conditions  and  the  expected  target 
action.  Some  of  the  motions  emulated  in  reference  6  include  fleeing  target,  random  walk,  front 
motion,  intercept,  shadow,  and  evasive  maneuver.  In  the  grid-based  method  of  reference  5,  the 
likelihood  of  signal  detection  (which  is  zero  for  target  positions  over  land)  automatically  zeros 
the  target  probability  density  in  the  constraint  region.  The  velocity  distribution  is  also 
maintained  on  a  grid,  and  its  evolution  is  driven  by  the  measurement  and  detection  likelihood. 


A  NON-HOMOGENEOUS  PROCESS  NOISE  MODEL 
FOR  TRACKING  NEAR  LANDMASS  CONSTRAINTS 

In  applying  modem  control  and  estimation  theory  to  the  tracking  of  ocean  targets,  the 
following  constant-velocity  motion  model  is  generally  applied: 

Xk+ 1  =  FkXk  +  Gk®k  >  (1) 

zk=hk(xk)  +  rik,  (2) 


where  xk  is  the  vector  representing  position  and  velocity  in  Cartesian  coordinates  at  time  4,  and 
zk  is  the  corresponding  measurement  vector.  Here,  the  plant  matrix  Fk  defines  the  deterministic 
component  of  the  target  state  transition  for  time  index  k  to  k  +  1 ,  and  it  is  solely  a  function  of  the 
time  increment  (4+1  -  4).  The  process  noise  term  cok  accounts  for  random  errors  and  stochastic 
disturbances  (maneuvers)  in  the  target  motion  description,  and  Gk  is  usually  considered  a 
constant  matrix.  The  measurements  zk  can  be  nonlinearly  related  to  the  target  state  (e.g.,  bearing 
and  range),  although  in  the  simulations  that  follow  in  this  report  direct  observations  of  target 
position  in  Cartesian  coordinates  were  used.  Both  the  measurement  noise  77  and  process  noise  co 
are  typically  considered  to  be  zero-mean  Gaussian  random  processes  with  fixed  covariance 
matrices  R  and  Q,  respectively. 

Classic  applications  of  Kalman  filter  theory  for  target  tracking  employ  homogeneous, 
stationary,  Gaussian,  white  noise  sequences.  In  contrast,  the  method  to  be  described  here  uses 
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non-homogeneous  process  noise  and  is  currently  being  investigated  at  the  Naval  Undersea 
Warfare  Center  (NUWC)  Division,  Newport,  RI,  for  use  when  tracking  in  the  presence  of 
landmass  constraints.  Traditionally,  a  mean  noise  value  is  lumped  into  the  deterministic  model, 
either  as  a  control  input  or  in  the  state  transition  function.  In  the  development  that  follows,  it 
will  be  instructive  to  consider  process  noise  models  with  nonzero  mean. 

It  is  assumed  here  that  the  waterborne  targets  of  interest  are  not  able  to  traverse  land 
barriers.  In  addition,  it  is  assumed  that  when  these  targets  are  in  the  vicinity  of  land  they  have  a 
greater  tendency  to  maneuver  (to  avoid  running  aground).  The  first  assumption  can  be 
approximately  modeled  in  the  target  motion  equations  by  a  large  coefficient  of  friction  over  the 
landmass  regions.  The  effects  of  a  frictional  force  can  be  included  in  the  state  transition  function 
or,  as  will  be  seen,  as  a  mean  on  the  process  noise  that  is  parametrically  described  in  terms  of  the 
current  target  velocity.  The  tendency  for  the  target  to  maneuver  can  be  modeled  by  an  elevated 
value  of  the  process  noise  variance.  As  a  result,  the  desired  target  motion  model  in  the  vicinity 
of  land  can  be  approximated  using  state  dependent  (non-homogeneous)  process  noise.  The 
problem  remains,  then,  to  characterize  this  state  dependence  in  a  manner  that  provides  a  smooth 
transition  from  open  ocean  to  near-land  conditions.  This  characterization  is  achieved  in  a  simple 
and  effective  manner  by  using  a  Gaussian  sum  representation  of  the  landmass  regions.  (Streit 
and  Luginbuhl  (reference  7)  applied  Gaussian  sum  representations  to  the  multitarget/multisensor 
data  association  problem,  and  their  work  inspired  such  application  to  the  present  problem.) 

Because  the  primary  reason  for  using  process  noise  in  this  instance  is  to  model  the  velocity 
bias  and  maneuver  uncertainty,  in  the  discussion  that  follows  the  process  noise  is  considered  to 
be  concentrated  in  velocity.  The  impact  on  position  is  achieved  through  the  integration  implicit 
in  the  state  transition  function.  For  the  application  at  hand,  one  can  consider  the  process  noise  to 
arise  from  multiple  sources,  which  will  be  termed  maneuver  influence  sources.  The  relative 
contribution  of  each  of  these  sources  depends  on  the  target  location.  One  maneuver  influence 
source  will  be  defined  for  the  open  ocean  and,  since  it  will  be  small  relative  to  the  other  sources 
and  negligible  in  their  presence,  it  will  be  considered  to  have  constant  weight  throughout  the 
region  of  interest.  In  addition  to  the  open-ocean  source,  there  is  a  maneuver  influence  source 
associated  with  each  of  the  components  in  a  mixture  representation  of  land.  These  sources 
receive  a  weight  that  is  proportional  to  a  target’s  proximity  to  the  particular  landmass  component 
in  question;  i.e.,  these  sources  are  weighted  by  an  exponential  form  directly  analogous  to  a 
Gaussian  probability  density  with  mean  and  covariance  determined  by  that  particular  component 
of  the  landmass  representation.  Finally,  these  weights  are  normalized  so  that  the  sum  over  all 
sources  is  unity;  see  equation  (3).  For  simplicity  and  tractability,  each  of  these  process  noise 
sources  is  considered  to  be  Gaussian,  with  specified  mean  and  covariance.  The  open-ocean 
source  is  modeled  as  zero  mean  (p.0  =  0)  and  with  a  relatively  small  covariance  typical  of  that 
used  in  classic  application  of  Kalman  filter  equations  to  the  constant- velocity  target  tracking 
problem.  The  sources  associated  with  landmass  components  are  modeled  with  a  mean  value 
emulating  a  frictional  force  that  prevents  motion  on  land  and  a  small  reflective  velocity  tailored 
to  model  any  land  avoidance  tendency.  The  covariance  of  these  land  sources  is  relatively  high, 
indicative  of  a  strong  tendency  to  maneuver.  Thus,  for  land  mode  /(  /  =  1,2, ...,  n),  the 
covariance  is  fixed  at  0,  and  the  mean  is  given  by 
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(3) 


Hi(v)  =  -av+vi, 

where  v  is  the  target  velocity  vector,  the  form  of  v,  is  designed  to  model  whatever  land- 
avoidance  action  is  expected  of  the  target,  and  a  corresponds  to  the  coefficient  of  friction.  In  the 
simulation  results  that  follow,  i>,  is  simply  a  small  velocity  vector  directed  away  from  the  land 
mode  in  question  and  toward  the  open  ocean.  It  is  important  that  the  coefficient  a  be  close 
enough  to  unity  that  a  target  approaching  land  slows  sufficiently  in  relatively  few  update 
intervals.  In  tests  performed  to  date,  the  value  a  =  1  has  been  used. 

The  weighted  sum  of  Gaussian  random  variables  is  also  Gaussian  distributed  with  a  mean 
value  given  by  the  weighted  sum  of  the  individual  means  and  a  covariance  given  by  the  weighted 
sum  of  the  individual  covariances. 


Putting  all  of  the  foregoing  elements  together  gives  the  following  functions  for  the  mean 
and  covariance  of  the  Gaussian  process  noise  model: 

0>k  =  N[n(xlk),Q(xtk)],  (4) 

fl(x)=tci(r)ni(v),  (5) 

/=0 

Q(x)  =  tc?(r)Qi,  (6) 

i=0 


where  xtk  =  (rlk,vlk)  is  the  estimated  target  state  (position  and  velocity)  at  time  4,  the  /t,are 
defined  in  equation  (3),  and  coefficients  C/r)  are  defined  as  follows: 


Q(r)  = 


^ 

7r0  +  t/r,  N(mj^j)\r 
;=I 


C,(r)  = 


jr,  N(mi,I.i)\r 
K0  +  t^jN(mj,I.j)\r 


i  =  1,  2,  ...,  n. 


(7a) 


(7b) 


Here,  N(m,,'Ld\r  signifies  the  multivariate  Gaussian  density  function,  with  mean  W;  and 
covariance  matrix  I,  evaluated  at  the  point  r.  It  is  sufficient,  for  the  purposes  of  this  formulation, 
to  equate  Kj  to  unity  for  i  >  0,  i.e.,  K\  =  7T2  =  .. .  n„  =  1 .  The  sum  appearing  in  the  denominator  of 
equation  (7), 

m=ixjN(mj,I.j)  \r,  (8) 

;'= i 
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corresponds  to  the  mixture  representation  of  the  landmass  region  (see  figure  1).  The  contour 
parametrically  defined  by  /(r)  =y  is  designed  to  approximate  the  shoreline  in  the  region  of 
interest  (figure  lb);  the  values  of  (mh  I„  i  =  1, 2, n)  and  n  are  design  parameters  to  be 
specified  for  each  region  of  interest.  The  values  of  7t0  and  y  are  design  parameters  that  are 
expected  to  be  fixed  for  all  regions.  An  illustration  of  the  process  noise  field  as  a  function  of 
location  is  shown  in  figure  2  for  the  land  boundary  shown  in  figure  1  and  a  hypothesized  target 
moving  in  the  positive  x-direction. 

The  descriptive  parameters  for  figure  2  are  given  in  table  1 .  The  figure  shows  a  grid 
sampling  of  the  process  noise  mean  vector  field  as  a  function  of  location;  centered  at  the  base  of 
each  vector  is  the  ellipse  representing  the  corresponding  process  noise  covariance  matrix.  Over 
navigable  water,  the  process  noise  mean  and  variance  are  approximately  zero;  they  become 
significant  only  near  land,  where  the  maneuver  influence  sources  are  appreciable. 


Figure  2.  Samples  of  the  Process  Noise  Field  Showing  Mean  Vectors  and 
Covariance  Ellipses  as  a  Function  of  Location 


Table  1.  Simulation  Parameters  for  Land  Parameters  and  Process  Parameters 


Although  the  target  kinematic  model  has  been  augmented  to  account  for  the  presence  of  a 
landmass,  the  target  density  promoted  by  the  Kalman  tracking  filter  is  still  Gaussian  in  nature. 
The  simplest  method  for  including  the  land  constraint  in  the  target  probability  density  is  to 
employ  a  window  function  that  is  unity  over  water  and  zero  over  land.  The  function  C0(r) 
(equation  7a)  is  a  smooth  function  that  closely  satisfies  the  above  condition  (see  figure  3). 

Thus,  a  simple  approximation  of  the  desired  target  density  is  provided  by 

p(x/z)  =  k  C0(r)  N(x,P)  \x,  (9) 

where  xT  =  (rT,vT),  A:  is  a  normalization  constant  such  that  p{xJ z)  integrates  to  1,  and  x,  P  are 
the  Kalman  filter  estimates  of  the  target  density  mean  and  covariance,  respectively.  Note  that  if 
constraint  regions  are  absent  in  the  region  of  interest,  this  representation  reverts  to  a  simple 
Gaussian  representation  of  the  target  density. 
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Figure  3.  Perspective  Plot  of  the  Window  Function  C  fr) 


SIMULATION  RESULTS 


The  formulation  presented  in  the  preceding  section  provides  a  compact  representation  of 
the  constraint  region  in  terms  of  a  sum  of  Gaussian  distributions,  a  smoothed  window  function 
for  the  target  density,  and  a  simple  and  effective  means  for  implementing  a  non-homogeneous 
process  noise  model.  Associated  with  each  term  of  the  Gaussian  sum  representation  of  land  are 
corresponding  process  noise  parameters  T]  and  O.  These  parameters  are  most  readily  stored  as 
velocity  components  of  the  land  element  state;  i.e., 


0,.  ={m, 


Q.= 


X  o 
o  a 


(10) 


Thus,  each  land  element  has  a  data  structure  not  unlike  that  of  the  target,  with  mean  vector  and 
covariance  matrix;  however,  the  information  is  used  in  a  much  different  manner. 

Several  simulations  were  performed  to  illustrate  the  algorithm’s  performance  under  varying 
data  conditions.  Coordinates  were  normalized  for  convenient  land  mixture  and  target 
representation;  time  was  similarly  scaled  to  keep  the  simulation  compact.  A  “data-rich”  scenario 
was  used  to  illustrate  tracking  through  land-avoidance  maneuvers;  a  “sparse-data”  scenario  was 
used  to  illustrate  the  track  maintenance  function;  a  “no-data”  (loss  of  contact)  scenario  was  used 
to  illustrate  the  ability  to  predict  the  probability  density.  In  all  cases,  the  land  model  delineated 
by  table  1  was  used. 
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Figure  4  illustrates  the  simulation  under  data-rich  conditions,  where  noisy  measurements  of 
position  are  available  at  each  time  index;  i.e., 

Z-k  =  H/s Xk  Vk’  (1 1) 

where  rjk  is  the  zero-mean  Gaussian  measurement  noise  at  time  tk,  and  rjk  ~  NfO.RJ.  In  figure  4a 
the  noisy  measurements  and  simulated  target  track  are  superimposed  on  the  plot  of  the  36  land 
mode  mean  positions  and  the  resulting  shoreline  contour.  The  parameters  pertaining  to  the  target 
and  measurement  process  are  given  in  table  2.  In  all  of  the  simulations  presented  here,  the 
algorithm  was  initialized  with  the  true  target  state  at  time  zero  (x0)  and  a  moderate  covariance 
(Pq),  since  it  was  the  behavior  near  the  land  boundary  that  was  of  primary  interest.  Figure  4b 
illustrates  the  algorithm’s  output  (dotted  curve)  when  the  presence  of  land  is  not  accounted  for. 
Solution  uncertainty  contours  are  shown  for  the  initial,  middle,  and  final  time  solutions.  Figure 
4c  shows  the  algorithm’s  output  when  the  non-homogeneous  process  noise  model  of  this  report 
is  applied.  The  increase  in  the  algorithm’s  gain  resulting  from  the  land-induced  process  noise 
allows  it  to  better  track  through  the  land-avoidance  maneuver,  yielding  a  more  accurate  solution 
and  a  more  consistent  solution  uncertainty  estimate. 


Figure  4a.  Data-Rich  Scenario 
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Figure  4b.  Algorithm  Output  When  No  Land  Constraint  Is  Applied 
for  Data-Rich  Scenario 


Figure  4c.  Algorithm  Output  Using  N on-Homogeneous  Process  Noise  Model 

for  Data-Rich  Scenario 


Table  2.  Simulation  Parameters  for  Target  Process  and  Measurement  Process 


Target  Process: 

rQ=  (-10, 3) 

(Initial  target  position) 

Vo7- =(1,0) 

(Initial  target  velocity) 

v„r=(0,l) 

(Instantaneous  maneuver  at  r  =  11) 

xj=  (r0T,  vf) 

(Position  and  velocity  state  vector) 

dt-  1.0 

Po  = 


0  0 

0  0 

0  0  0.01  0 
0  0.01 


1  0 
0  1 


0  0 


(Initial  target  covariance) 


'  1 

0 

dt 

O' 

dt 

o' 

Fk  = 

0 

0 

1 

0 

0 

1 

dt 

0 

Gk~ 

0 

1 

dt 

0 

_0 

0 

0 

1 

_  0 

1  _ 

Measurement  Process: 


Hk  = 


'1  0  0  0' 
0  10  0 


Rk 


1  O' 
0  1 


Figure  5  illustrates  the  same  scenario  under  sparse-data  conditions,  where  noisy  position 
measurements  are  available  every  seventh  time  index,  i.e.,  at  times  1,8,  15,  and  22.  Figure  5a 
shows  the  simulated  target  track  and  measurement  set.  Figure  5b  shows  the  algorithm’s  output 
when  land  is  not  accounted  for,  and  one  can  see  an  obvious  violation  of  the  land  constraint.  In 
figure  5c  it  is  seen  once  again  that  the  land-induced  process  noise  increases  the  algorithm’s  gain 
to  allow  tracking  through  the  land-avoidance  maneuver,  and  it  provides  a  consistent  estimate  of 
the  solution  uncertainty.  The  consistency  of  the  uncertainty  estimate  is  highlighted  by  including 
the  uncertainty  contours  at  time  21,  just  prior  to  processing  the  last  position  measurement. 
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Figure  5c.  Algorithm  Output  Using  N on-Homogeneous  Process  Noise  Model 

for  Sparse-Data  Scenario 


The  predictions  for  a  lost  contact  (no  data)  are  illustrated  in  figure  6  where  the  same 
scenario  used  in  figures  4  and  5  is  simulated  without  any  available  measurements.  Because  of 
the  perfect  initialization,  the  simulated  and  predicted  tracks  of  figure  6a  are  coincident  for  the 
first  half  of  the  scenario.  However,  in  the  absence  of  a  land-constraint  model,  the  need  for  a 
land-avoidance  maneuver  is  not  accounted  for  in  the  last  half  of  the  predicted  target  track. 
Figure  6b  illustrates  how  the  non-homogeneous  process  noise  field  more  accurately  models  the 
actual  estimate  uncertainty  while  also  enforcing  the  land  constraint  on  target  motion.  Note  how 
the  final  time  contours  of  the  windowed  density  wrap  around  the  landmass.  Since  the  target 
could  have  chosen  any  land-avoidance  maneuver — and  none  was  more  probable  than  another 
(t>, «  0) — little  movement  of  the  mean  location  is  predicted  after  contact  with  shore.  The 
predicted  location  uncertainty  is,  however,  large  and  representative  of  the  lack  of  maneuver 
information. 
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Figure  6a.  Algorithm  Output  When  No  Land  Constraint  Is  Applied 

for  No-Data  Scenario 


Figure  6b.  Algorithm  Output  Using  N on-Homogeneous  Process  Noise  Model 

for  No-Data  Scenario 


CONCLUSIONS 


As  has  been  demonstrated  by  the  simulation  results  presented  in  this  report,  the  net  effect 
of  applying  an  involved  model  of  the  target  process  noise  field  to  the  target  tracking  algorithm  is 
that  new  data  are  readily  incorporated  in  the  manner  of  a  Kalman  filter,  and  the  altered  target 
behavior  as  it  approaches  a  constraint  boundary  is  accommodated  in  a  compact  form.  Properly 
used,  such  a  model  will  prevent  constraint  violation  and  will  allow  a  more  accurate  portrayal  of 
the  state  uncertainty  near  boundaries.  The  regions  of  constraint  violation  are  represented  using  a 
mixture  model  that  is  fixed  for  a  given  region  of  interest,  thus  limiting  the  burden  of  excessive 
database  interaction  during  filter  iterations. 

The  simulation  results  illustrate  the  shift  in  focus  from  predictions  of  the  mean  estimate 
when  in  open  water  to  predictions  of  the  estimate  uncertainty  growth  when  near  land.  This 
refocus  is  a  direct  result  of  the  assumption  that  an  intelligent  agent  is  directing  the  target  platform 
and  so  it  will  maneuver  to  avoid  running  aground.  The  lack  of  specific  maneuver  information  is 
manifested  in  the  increase  in  the  target  estimate  uncertainty.  In  the  tracking  model  presented 
here,  the  lack  of  maneuver  information  is  accounted  for  by  using  maneuver  influence  sources  in 
the  process  noise  formulation.  The  mode  shift  from  the  point  estimate  to  the  uncertainty  region 
is  an  important  aspect  of  tracking  near  constraint  boundaries. 

Although  quantitative  performance  measures  would  be  premature  at  this  juncture,  factors 
affecting  performance  include  the  number  of  components  representing  the  constraint  region,  the 
parameterization  of  those  components,  and  the  rate  of  algorithm  update.  Greater  accuracy  in  the 
a  posteriori  probability  estimate  may  be  achieved  by  using  multiple  components  for  the  target 
density  itself  (reference  8)  or  by  representing  it  by  a  multivariate  Gram-Charlier  series  or 
Edgeworth  expansion  (references  9-11).  The  efficacy  of  these  representations  will  be  the  subject 
of  future  investigations. 


17 


REFERENCES 


1 .  “CAAM  User's  Guide  -  Preliminary  Draft,”  Code  6011,  Naval  Undersea  Warfare  Center 
Detachment,  New  London,  CT,  November  1992. 

2.  “ASWTDA  Methodology  Guide,”  Sonalysts  Document  CA1-0501-029,  prepared  for  Naval 
Undersea  Warfare  Center  Detachment,  New  London,  CT,  17  July  1992. 

3.  S.W.  Tomasi,  P.  Kelly,  and  G.  Morrison,  “Functional  Requirements  Definition  for  Visual 
Integration  of  Scene  Information  Over  Networks  (VISION),  Preliminary  Draft,”  Naval 
Undersea  Warfare  Center  Division,  Newport,  RI,  22  May  1995. 

4.  Donald  A.  Pierre,  Optimization  Theory  With  Applications,  John  Wiley  and  Sons,  New  York, 
1969. 

5.  Lawrence  D.  Stone  and  Thomas  M.  Kratzke,  “Comparison  of  Linear  and  Non-Linear 
Trackers,”  Metron  Inc.  report  to  the  Naval  Research  Laboratory,  Washington,  DC,  16 
December  1991. 

6.  “Module  Design  Document  for  the  Combined  DEEP/Shallow  Water  Data  Fusion 
Algorithm,”  Sonalysts  Document  prepared  for  B.  I.  Incze  (Code  2111),  Naval  Undersea 
Warfare  Center  Detachment,  New  London,  CT,  7  May  1993. 

7.  Roy  L.  Streit,  and  Tod  E.  Luginbuhl,  “Probabilistic  Multi-Hypothesis  Tracking,” 
NUWC-NPT  Technical  Report  10,428,  Naval  Undersea  Warfare  Center  Division,  Newport, 
RI,  15  February  1995. 

8.  D.  L.  Alspach  and  H.  W.  Sorenson,  “Approximation  of  Density  Functions  by  a  Sum  of 
Gaussians  for  Nonlinear  Bayesian  Estimation,  ”  Proceedings  of  the  Symposium  on  Nonlinear 
Estimation  Theory  and  its  Applications ,  San  Diego,  CA,  21-23  September  1970. 

9.  Vincent  J.  Aidala,  "The  Multivariate  Gram-Charlier  Series  Applied  to  Random  Signal 
Detection,"  NUWC-NPT  Technical  Report  10,486,  Naval  Undersea  Warfare  Center 
Division,  Newport,  RI,  14  July  1995. 

10.  H.  W.  Sorenson  and  A.  R.  Stubberud,  “Non-linear  Filtering  by  Approximation  of  the  A- 
posteriori  Density,”  International  Journal  of  Control,  vol.  8,  no.  1,  1968,  pp.  33-51. 

11.  Peter  S.  Maybeck,  Stochastic  Models,  Estimation,  and  Control,  Volume  2,  Academic  Press, 
New  York,  1982. 


18 


INITIAL  DISTRIBUTION  LIST 


Addressee  No.  of  Copies 

Commander,  Submarine  Force  Atlantic  1 

Commander,  Submarine  Force  Pacific  1 

Commander,  Submarine  Development  Squadron  12  1 

Defense  Technical  Information  Center  12 

Advanced  Research  Projects  Agency  1 

Chief  of  Naval  Research  (ONR-333  (J.  Fein),  ONR-342  (T.  McMullen))  2 

Center  for  Naval  Analyses  1 

Naval  Postgraduate  School  1 

Naval  Sea  Systems  Command  (ASTO,  ASTO-E  (P.  Kelly),  PMO-41 1,  PMO-415)  4 

Naval  Surface  Warfare  Center,  Dahlgren  Division  (Code  N51  (M.  Williams))  1 


